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ABSTRACT 


In  an  attempt  to  find  optimum  reference  waveforms  for 
matched  filtering  long  period  signals,  several  reference  waveforms  varying 
in  duration  and  signal  contents  were  taken  from  Alaskan  Long  Period  Array 
recordings  of  major  events  in  the  Sinkiang  Province  area.  Each  of  these 
reference  waveforms  was  used  as  a  matched  filter  in  processing  a  number 
of  other  events  in  the  same  area.  The  events  differ  in  magnitude,  location 
and  depth.  Analysis  of  the  signal -to -noise  ratio  gains  (relative  to  bandpass 
filtering)  showed  that  optimum  processing  of  all  events  required  a  large 
number  of  reference  waveforms.  This  processing  yielded  up  to  8  dB  improve  - 
ment  over  to-date  routine  matched  filter  processing  of  Central  Asia  evenls. 
Several  other  aspects  of  reference  waveform  matched  filtering  also  are  discussed. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and 
this  document  is  subject  to  later  revision  as  ma;  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should  not  be  inter¬ 
preted  as  necessarily  representing  the  official  policies,  either  expressed  or 
implied,  of  the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical 
Applications  Center,  or  the  US  Government. 
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SECTION  I 


INTRODUCTION 


Previous  studies  (Texas  Instruments,  1971;  1972)  have  shown 
that  when  using  a  long  period  signal  from  one  seismic  event  as  a  reference 
waveform  for  matched  filtering  other  events  in  a  given  area,  some  uncer¬ 
tainties  exist  about  the  selection  of: 

•  the  reference  event 

•  the  part  of  the  reference  event  signal  that  will  yield  maximum 
processing  gain. 

This  study  analyzes  the  effects  that  various  choices  of  reference  waveforms  have 

/ 

on  matched  filter  processing  using  existing  routine  processing  programs.  Ana¬ 
lysis  is  performed  on  Alaskan  Long  Period  Array  (ALPA)  recordings  of  events  in 
the  Sinkiang  Province  area. 

Section  II  of  this  report  discusses  general  signal  processing  and 
seismological  considerations  regarding  reference  waveform  matched  filtering. 
Section  III  describes  the  analysis  approach  taken  on  the  study,  and  the  data  pro¬ 
cessing  procedures  followed.  In  Section  IV  the  results  are  interpreted  in  terms 
of  the  optimum  reference  event  and  the  optimum  reference  waveform  length, 
taking  into  account  the  effects  of  location  and  magnitude.  This  section  also  com¬ 
pares  optimum  processing  gains  with  to-date  routine  processing  gains,  and  ana¬ 
lyzes  gain  differences  between  the  vertical  and  the  radial  component.  Finally, 
Section  V  summarizes  the  results  of  the  study. 


SECTION  II 


GENERAL  CONSIDERATIONS 

A.  INTRODUCTION 

Before  starting  the  experimental  search  for  an  optimum 
type  reference  waveform  we  review  some  signal-processing  and  seismo- 
logiaal  elements  involved  in  reference  waveform  matched  filtering,  which 
form  the  basis  for  the  study  organization  and  the  data  analysis.  Also,  the 
matched  filtering  effects  on  false  alarm  rate  and  surface  wave  magnitude 
measurement  are  discussed  briefly. 

B.  SIGNAL  PROCESSING  CONSIDERATIONS 

For  the  case  of  a  noise-free,  uniform-envelope,  dispersed 
reference  waveform  which  is  an  exact  replica  of  the  target  signal  buried  in 
band-limited  white  background  noise,  the  processing  gain  approximately  equals 
the  time-bandwidth  product  (Capon,  1967).  When  using  the  signal  from  one 
seismic  event  to  process  another  event  in  the  same  area,  the  reference  wave¬ 
form  is  not  noise-free,  does  rot  have  a  uniform  envelope,  is  not  an  exact 
replica  of  th*e  other  signal,  and  the  noise  is  non-white.  Therefore,  the  ex¬ 
pected  processing  gain  in  this  type  of  filtering  is  less  than  the  time-bandwidth 
product.  However,  a  long-duration,  wide-band,  high  signal -to -noise  ratio  re¬ 
ference  waveform  still  would  yield  maximum  processing  gain  provided  that 
significant  signal  similarity  is  mair  lained  over  the  entire  signal  duration. 

When  signal  similarities  are  of  a  limited  duration,  a  shorter  duration  re¬ 
ference  waveform  must  be  selected.  Because  of  dispersion  this  means  that 
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the  reference  waveform  contains  a  narrower  frequency  band,  and  the  pro¬ 
cessing  gain  is  reduced  by  both  factors  in  the  time-bandwidth  product. 

Another  aspect  of  reference  waveform  processing  is  the 
effect  on  noise.  Considering  first  the  bandpass  effect  of  the  matched  filter 
we  note  that  the  shape  of  the  reference  waveform  magnitude  spectrum  (which 
generally  is  non-white)  can  change  the  output  noise  power  relative  to  the  input 
noise  power.  The  effect  can  either  increase  or  decrease  the  processing  gain, 
depending  on  the  relative  shapes  of  the  signal  and  noise  spectra.  Second’y,  if 
some  of  the  noise  waveforms  are  in  phase  with  part  of  the  reference  wave¬ 
form  the  noise  will  correlate  to  some  extent  with  the  reference  waveform, 
which  reduces  the  processing  gain. 

C.  SEISMOLOGICAL  CONSIDERATIONS 

Signal-to -noise  ratio  and  signal  similarities  of  events  are 

determined  by: 

•  Source  mechanism  (including  depth  effects) 

•  Event  magnitude 

•  Travel  path 

The  source  mechanism  determines  the  frequency  spectrum  and  the  radiation 
pattern.  Depending  on  the  geological  structure,  events  at  different  locations 
and  at  different  depths  in  general  have  different  source  mechanisms.  Events  of 
nearly  the  same  location  generally  have  similar  mechanisms,  but  the  event 
magnitudes  may  differ.  Signals  generated  at  different  locations  propagate 
along  different  paths,  causing  differences  in  dispersion,  multipathing,  and 
frequency-dependent  attenuation.  The  spectrum  generated  at  the  source,  the 
dispersion,  the  multipathing,  and  the  frequency-dependent  attenuation  deter¬ 
mine  the  signal  shape  upon  arrival  at  thearray;  the  magnitude,  the  radiation 
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pattern,  the  attenuation,  and  the  noise  behavior  during  the  signal  arrivals 
determine  its  signal -to -noise  ratio. 

In  reference  waveform  matched  filtering  the  reference 
events  are  chosen  for  their  expected  signal  similarity  and  signal-to-noise 
ratio  based  on  event  location,  magnitude  and  depth. 

D.  FALSE  ALARM  CONSIDERATIONS 

Since  a  relatively  large  number  of  reference  waveforms 
must  be  generated  to  enable  finding  optimum  waveforms,  this  could  increase 
the  false  alarm  rate,  depending  on  the  detection  criterion  used.  In  ALPA 
and  NORSAR  long  period  processing  the  criterion  requires  that,  within  the 
expected  signal  arrival  gate,  the  matched  filter  output  contain  a  peak  at  least 
3  dB  larger  than  any  noise  peak  within  about  30  minutes  of  the  signal  gate 
noise.  A  false  alarm  is  created  when,  according  to  this  criterion,  a  signal 
is  detected  in  the  filter  output  while  the  input  consists  of  noise  only. 

The  f  alse  alarm  rate  can  be  estimated  by  applying  matched 
filter  reference  waveforms  to  a  number  of  noise  gates.  The  false  alarm  rate, 
when  using  the  above  detection  criterion,  is  expected  to  be  less  than  1  percent, 
even  if  more  than  one  reference  or  chirp  waveform  are  used  (Texas  Instru¬ 
ments,  1973).  For  this  study  the  false  alarm  rate  maybe  somewhat  higher 
because  of  the  large  number  of  reference  waveform  used.  Time  did  not 
permit  an  analysis  of  the  false  alarm  rate;  this  should  be  done  in  the  future. 
We  would  not  expect  it  to  be  significantly  increased,  however. 

E.  SURFACE  WAVE  MAGNITUDE  CONSIDERATIONS 

The  surface  wave  magnitude  is  given  by 

M  =  log  +  1. 66  log  /\ 

S  1 
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(ii.  i) 


where  A  =  signal  peak  amplitude  (m n  ) 

T  =  signal  peak  period  (sec) 

A  =  epicentral  distance  (degrees) 

Matched  filtering  generally  will  cause  changes  in  A  and  T  which  would  distort 
the  Mg  measurement.  In  ALPA  and  NORSAR  long  period  processing  the 
filter  gain  is  compensated  by  subtracting  the  logarithm  of  the  regional  average 
gain  a  s  follows: 

A  1 

M  =  log  —  +  1 .  66  log  A  -  — -  G  (II.  2) 

s  T  20  m 

where  G^  =  regional  average  filter  gain  in  dB. 

A  gain  deviation  of,  for  example,  3  dB  about  the  average 

gain  would  cause  a  deviation  in  Mg  of  0.15.  A  change  in  signal  peak  period 

T  of,  for  example,  5  sec  would  cause  an  M  deviation  of  0.  10.  Deviations 

s 

in  A  and  T  are  in  general  smaller  and  the  resulting  change  of  less  than  0.  25 
in  Mg,  compared  with  a  typical  standard  deviation  of  0.  5,  is  not  considered 
a  major  distortion  in  the  measurement  of  surface  wave  magnitudes. 
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SECTION  III 


ORGANIZATION  AND  PROCEDURES 

A.  INTRODUCTION 

To  pursue  the  objective  of  this  study,  three  Sinkiang  Province 
events  were  chosen  as  reference  events.  From  each  reference  event  signal 
trace,  eight  different  reference  waveform  lengths  were  generated  by  selecting 
different  time  windows.  Each  of  the  total  of  twenty-four  reference  waveforms 
thus  obtained  was  then  used  as  a  matched  filter  in  processing  each  of  twenty- 
one  Sinkiang  Province  test  events.  The  resulting  signal -to-noise  ratio  gains 
relative  to  bandpass  filtering  were  compared  and  analyzed. 

The  Sinkiang  Province  area  consists  of  mountains,  rivers  and 
desert.  The  majority  of  the  twenty-one  events  were  located  in  the  mountainous 
region  on  the  Soviet-China  border;  the  location  of  the  events  is  shown  in  Fi¬ 
gure  III  - 1  by  code  number  as  listed  in  Table  III  - 1  which  groups  the  events  ac¬ 
cording  to  their  distance  to  the  three  reference  events.  The  table  lists  the 
event  code  number,  the  group  number,  the  event  name,  the  event  latitude  and 
longitude,  the  great  circle  distance  and  azimuth  from  array  to  event  epicenter, 
the  source  depth,  the  body-wave  magnitude  and  data  source,  the  number  of  ,ood 
sites  used  in  array  beaiusteering,  and  the  distance  between  test  event  and  each 
reference  event.  The  reference  events  (168,  170,  002)  are  underlined.  Data 
sources  were  given  by  previous  processirg  and  reports,  and  are  abbreviated 
as  follows: 

P  =  Preliminary  Determination  of  Epicenter  (P.  D.  E.) 

L  =  Large  Aperture  Seismic  Array  (L.  A.  S.  A.) 
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N  -  Norwegian  Seismic  Array  (NORSAR) 

I  =  International  Seismic  Bulletin  (I.  S.  B.) 

M  -  International  Seismological  Month  (I.  S.  M. ) 
Blank  =  non-P.  D.  E.  ,  but  otherwise  not  known 


The  test  events  of  group  I  are  within  61  km  of  the  reference 
events  refs.  168,  170;  their  distance  to  ref.  002  is  about  450  km.  The  events 
in  group  II  are  closest  to  ref.  002;  their  distance  to  refs.  168,  170  varies  from 
236  to  684  km.  The  events  of  groups  III  and  IV  are  all  located  closer  to  refs. 
168,  170  than  to  ref.  002,  but  their  distance  to  refs.  168,  170  is  considerab-y 
greater  than  those  of  the  group  I  events.  For  group  III  these  distances  are 
117  to  279  km,  for  group  IV  the  distances  vary  from  446  to  523  km.  Group  V 

is  comprised  by  only  one  event  having  approximately  equal  distance  to  each  of 
the  three  reference  events. 


B. 


REFERENCE  EVENTS 


To  meet  tins  study's  objective  we  need  a  set  of  reference  events 
with  different  locations  and  magnitudes,  and  with  a  sufficient  signal -to-noise 
ratio.  The  following  set  of  reference  events  satisfied  these  requirements: 


ref.  002: 
ref.  168: 
ref.  170: 


SIN-002-10AL,  m  =  5.  2 

b 

SIN/163/15AL,  m  =  4.  9 

b 

SIN/170/17AL,  m  =  5.2 

b 


Refs.  168  and  170  were  located  24  km  from  each  other  and  within  6l  km  from 
any  event  in  group  I  (see  Table  III-l).  Ref.  002  was  located  430  km  from  refs. 
168,  170  and  within  270  km  frcm  any  event  in  group  II.  The  depths  of  refs.  168 
170  were  33  km,  ref.  002  had  19  km  depth.  Because  of  the  location  and  depth 
differences  between  refs.  168,  170  on  the  one  hand  and  ref.  002  on  the  other, 
r ef'  002  most  likelY  had  a  source  mechanism  different  from  that  of  refs.  168, 
170.  The  reference  event  bandpass  signal  traces  are  shown  in  the  Figures  III-: 
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III-3,  III -4,  which  also  illustrate  the  reference  waveform  lengths  selected  as 
described  in  the  next  subsection.  The  signals  from  ref.  163  and  170  are  indeed 
similar  in  shape  (although  different  in  signal -to -noise  ratio),  while  the  ref.  002 
signal  is  different  in  both,  signal  shape  and  signal-to -noise  ratio. 

C.  REFERENCE  WAVEFORMS 

In  order  to  search  for  the  part  of  a  reference  event  signal 
trace  that  will  give  maximum  processing  gain  when  used  as  a  matched  filter 
reference  waveform,  eight  different,  but  overlapping,  parts  were  selected 
from  each  reference  event  signal  trace,  for  each  of  the  beamsteered  compo¬ 
nents.  Despite  the  signal  differences  between  the  three  reference  events,  the 
reference  waveforms  can  be  roughly  characterized  as  follows,  see  Figures 
III-2,  III- 3,  III  -  4: 

length  A:  first  significant  wa'  etrain  after  arrival  time 

length  B:  A  plus  second  wavetrain 

length  C:  B  plus  third  wavetrain 

length  D:  C  plus  fourth  wavetrain 

length  E:  D  plus  remainder  of  signal  arrival 

length  X:  wavetrain  containing  maximum  signal 

length  Y:  X  plus  one  adjacent  wavetrain 

length  Z:  X  plus  two  adjacent  wavetrains 

A  wavetrain  in  this  instance  is  considered  a  smoothly  modulated  signal  portion 
between  "nodes"  or  between  significant  phase  changes.  Because  of  the  difference 
in  signal  shape  of  the  three  reference  events  the  above  characterization  cannot 
be  very  rigid;  for  example  type  A  of  ref.  002  has  different  signal  content  and 
will  perform  differently  than  type  A  of  ref.  170. 

Because  of  dispersion,  length  A  in  general  contains  the  low 
frequency,  early  arrival  part  of  the  signal.  Lengths  B,  C,  D  and  E  include 
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increasingly  higher  frequencies  and  later  arrivals,  with  length  E  containing 
the  entire  signal,  including  multiple  arrivals.  Length  X  is  a  short  duration, 
center-frequency-band  waveform  :ontaining  the  highest  signal  power:  lengths 
Y  and  Z  are  extensions  of  the  X  waveform  . 

The  interval  specification  of  a  reference  waveform-J ,o ve 
wave  component  is  independent  of  the  Rayleigh-reference  waveform  speci¬ 
fication.  The  two  Rayleigh-reference  waveform  components  are  formed  by 
identical  interval  specification. 

D.  DATA  PROCESSING 

1.  Introduction 

Most  data  processing  was  done  with  existing  programs  used 
in  routine  processing  (Texas  Instruments,  1972).  These  programs  perform 
frequency  domain  rather  than  time  domain  filtering  because  greater  speed  is 
achieved  with  the  use  of  fast  Fourier  transforms.  The  various  processing 
steps  are  described  briefly  below. 

2.  Event  Signals 

All  signals  used  in  this  study  were  previously  edited  for 
spike  removal  and  for  elimination  of  bad  sites.  The  resulting  signals  were 
rotated  to  yield  the  transverse  Love  wave  component  (LQ-T),  and  the  ver¬ 
tical  and  radial  Rayleigh  wave  components  (LR-V,  LR-R).  The  array  signals 
were  obtained  by  beamsteering  all  available  sites  (Texas  Instruments,  1972). 
Before  matched  filtering,  the  test  event  signals  were  filtered  through  the 
0.025-0.  55  Hz  frequency  band  (40  to  18  seconds),  the  reference  waveforms 
were  not  bandpassed. 
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3. 


Reference  Waveforr  *  Matched  Filtering 


Reference  waveform  lengths  were  specified  as  described 
previously.  Where  necessary,  the  waveform  was  extended  with  zero  values 
to  obtain  a  box  car  window  with  a  total  of  2048  2 -sec  samples.  This  signal 
was  Fourier  transformed,  and  the  complex  conjugate  taken.  The  use  of 
other  windows  w as  not  considered  in  this  study. 

Matched  filtering  was  performed  by  taking  the  Fourier  trans¬ 
form  of  the  test  event  signal  after  extending  it  with  zero  values  to  2048  2-sec 
samples,  and  multiplying  the  complex  spectrum  with  the  complex  conjugate 
of  the  reference  waveform  spectrum.  The  result,  the  complex  cross  power 
spectrum,  then  was  transformed  back  into  the  time  domain  to  yield  the  matched 
filter  output. 

4.  Signal-to-noise  Ratio  Computation 

For  long  period  array  evaluation  projects  the  true  signal-to- 
noise  ratio  (SNR)  is  defined  as  the  ratio  between  the  maximum  value  of  signal 
plus  noise  during  the  expected  signal  arrivals  and  the  RMS  noise  value  taken 
over  a  time  gate  just  prior  to  the  signal  arrival  (Texas  Instruments,  1972). 

The  SNR  was  computed  automatically,  however  output  traces  obtained  by 
CALCOMP  plotting  were  used  to  check  the  results. 

5.  Processing  Gain  Computation 

For  each  test  event  the  SNR  of  the  bandpass  filter  output  is 
computed  first,  then  the  SNR  is  determined  for  each  of  the  eight  matched  filter 
outputs  per  component,  lor  each  reference  event.  The  SNR  gain  (GN)  of 
matched  filtering  over  bandpass  filtering  is  obtained  as  the  ratio  of  matched 
filter  output  SNR  to  the  bandpass  filter  output  SNR.  Finally,  the  matched 
filter  processing  gain  over  bandpass  filtering  is  expressed  in  decibels: 

GDB  =  20  log  GN. 


SECTION  IV 


ANALYSIS 


A.  INTRODUCTION 

This  section  presents  and  analyzes  the  data  'rocessing 
results  for  optimum  reference  event  and  optimum  reference  waveform 
characteristics.  Detailed  numerical  results  are  presented  in  tabular  form 
for  each  test  event  in  Appendix  A.  The  m  ults  of  all  test  events  are  combined 
in  the  Figures  IV-1,  IV-2,  IV-3  and  in  Table  IV-1. 

Figures  IV-1,  IV-2,  IV-3  show  the  matched  filter  processing 
gains  per  component  over  bandpass  filtering,  obtained  by  the  reference  waveform 
lengths  A,  B,  C,  D,  E,  X,  Y,  Z,  in  that  order,  of  each  reference  event  as 
applied  in  processing  the  twenty -one  test  events.  Table  IV-1  presents  the 
optimum  gain  over  bandpass  filtering,  and  indicates  which  are  the  optimum 
or  near-optimum  reference  waveform  lengths  per  component,  for  each  reference 
event.  To  facilitate  the  analysis  the  most  relevant  parameters  from  Table 
II-l  are  repeated,  and  the  bandpass  signal -to -nois e  ratios  are  included  in 
the  table.  Gains  of  less  than  2  dB  are  omitted  in  the  optimization  analysis. 

In  the  analysis  we  will  first  search  the  results  for  the  reference 
events  that  are  optimum  for  some  set  of  event-.  From  each  optimum  reference 
event  we  then  will  try  to  find  the  optimum  reference  waveform  characteristics. 
Finally,  the  relationships  between  the  processing  gain  and  the  parameters 
distance,  magnitude  and  signal-to-noise  ratio  are  investigated,  the  optimum 
gams  are  compared  with  routine  processing  gains,  and  the  matched  filter  be¬ 
havior  of  the  two  Rayleigh  wave  components  is  discussed. 
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FIGURE  IV-l-b 

PROCESSING  GAINS  VERSUS  REFERENCE  WAVEFORM  LENGTHS, 

REF.  002,  LR-V 
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FIGURE  IV- 1-  c 

PROCESSING  GAINS  VERSUS  REFERENCE  WAVEFORM  LENGTHS, 

REF.  002,  LR-R 
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FIGURE  IV-2-b 

PROCESSING  GAINS  VERSUS  REFERENCE  WAVEFORM  LENGTHS, 

REF.  168,  LR-V 
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FIGURE  IV-2-c 

PROCESSING  GAINS  VERSUS  REFERENCE  WAVEFORM  LENGTHS 

REF.  168,  LR-R 
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FIGURE  IV- 3 -a 

PROCESSING  GAINS  VERSUS  REFERENCE  WAVEFORM  LENGTHS, 

REF.  170,  LQ-T 
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FIGURE  IV-3-b 

PROCESSING  GAINS  VERSUS  REFERENCE  WAVEFORM  LENGTHS, 

REF.  170,  LR-V 
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FIGURE  IV-3-c 

PROCESSING  GAINS  VERSUS  REFERENCE  WAVEFORM  LENGTHS, 

REF.  170,  LR-R 
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B. 


OPTIMUM  REFERENCE  EVENTS 


In  the  search  for  optimum  reference  events  we  scan  the 
results  in  Table  IV-1  and  Figures  IV-1,  IV-2,  IV-3.  Starting  with  the  rela¬ 
tively  small  area  containing  the  events  of  group  I,  we  see  that  for  most  events 
either  ref.  168  or  ref.  170  yields  maximum  processing  gain.  The  exceptions 
are  event  281*09  which  correlates  better  with  ref.  002  for  all  three  com¬ 
ponents,  and  event  170/21  for  which  the  transverse  component  matches  better 
with  a  ref.  002  waveform.  The  only  differences  between  event  281*09  and  the 
other  group  I  events  are  its  lower  magnitude  and  its  20  to  35  km  larger  dis¬ 
tance  to  refs.  168,  170;  it  may  also  have  a  different  depth.  The  transverse 
component  of  event  170/21  is  lacking  in  high  frequency  energy  and  so  its 
signal  does  not  match  well  with  the  refs.  168,  170  waveforms.  The  Rayleigh 
wave  components  of  this  event  correlate  well  with  these  reference  waveforms. 
Gain  differences  between  ref.  168  and  ref.  170  in  general  are  less  than  1.  5  dB. 

For  most  group  II  events,  all  located  closer  to  ref.  002  than 
to  refs.  168,  170,  the  optimum  event  is  ref.  002.  One  exception  is  the  low 
signal  -  to -noi  s  e  ratio  event  037-07  which  shows  an  overall  low  or  even  negative 
gain  compared  with  bandpass  filtering.  The  second  exception  is  the  vertical 
component  of  event  221*01,  showing  a  1.  2  dB  better  match  with  ref.  168 
waveforms.  (The  radial  component  is  1.4  dB  better  for  ref.  002,  and  so  the 
Rayleigh  wave  gains  are  comparable  for  the  two  reference  waveforms.) 

For  the  group  III  events,  located  closer  to  refs.  168,  170 
than  to  ref.  002,  either  ref.  168  or  ref.  170  constitutes  the  optimum  reference 
event,  with  less  than  1.  5  dB  gain  differences.  Only  the  transverse  component 
of  event  154-06  shows  a  slightly  better  correlation  (0.  8  dB)  with  ref.  002. 

The  matched  filtering  performance  of  the  events  in  groups 
IV  and  V  is  too  poor  to  be  considered  in  the  optimization  analysis;  these 
events  are  all  more  than  400  km  from  any  reference  event. 
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Thus,  we  may  conclude  that  distance  seems  the  most  impor¬ 
tant  factor  in  selecting  a  reierence  event,  and  in  general  it  should  be  chosen 
closest  to  the  events  to  be  processed.  Even  within  a  relatively  small  area, 
however,  an  event  may  occasionally  correlate  better  with  a  reference  event 
located  at  a  much  greater  distance.  Such  cases  seem  unpredictable.  Further¬ 
more,  reference  events  in  approximately  the  same  location,  but  of  different 
magnitude,  seem  to  perform  equivalently,  which  would  be  expected  provided 
the  reference  events  were  relatively  large. 

C.  OPTIMUM  REFERENCE  WAVEFORMS 

Once  an  optimum  reference  event  is  established  one  must 
determine  the  part  of  its  signal  trace  which  constitutes  a  reference  waveform 
that  yields  maximum  matched  filter  processing  gain.  If  for  each  optimum 
reference  event  the  optimum  waveform  length  is  similar,  we  have  established 
an  optimum  length  reference  waveform.  If  we  find  different  reference  wave¬ 
form  lengths,  we  may  be  able  to  establish  some  characteristics  which 
facilitate  the  selection  of  optimum  waveforms. 

To  pursue  the  above  objective  we  again  scan  Table  IV- 1  and 
Figures  IV-1,  IV-2,  IV-3,  and  focus  on  the  reference  waveform  types  yielding 
optimum  or  near-optimum  matched  filter  processing  gain.  Near-optimum 
gain  in  this  instance  means  a  gain  that  is  within  2  dB  of  the  maximum  gain. 
Since  the  transverse  component  reference  waveform  is  defined  independently 
from  the  two  Rayleigh  wave  component  reference  waveforms  their  optimum 

waveform  lengths  may  differ.  The  two  Rayleigh  wave  reference  waveforms, 
however,  must  be  considered  simultaneously,  since  they  are  defined  by  the 
same  time  interval  specification. 

In  the  previous  subsection  it  was  shown  that  for  the  events 
of  groups  I  and  III  combined,  with  the  exception  of  event  281*09  and  the 
transverse  component  of  event  170/21,  either  ref.  168  or  ref.  170  constitutes 
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the  optimum  reference  event.  Disregarding  these  exceptions  we  will  now 
investigate  if  for  either  reference  event  there  is  a  reference  waveform  length 
that  is  commonly  optimum  for  all  events  of  this  set.  Thus,  for  transverse 
component  processing,  we  consider  the  ten  events  167/13,  168/15,  170/17, 

184*04,  181*13,  042-05,  047-23,  064-04,  080-21,  154-06.  For  Rayleigh 
wave  processing  we  consider  the  same  events  plus  event  170/21,  a  total  of 
eleven  events.  When  considering  ref.  1 6 8  processing  vve  exclude  the  processing 
results  for  ref.  1 68  matched  to  itself.  Similarly,  we  omit  ref.  170  matched 
to  itself.  Thus,  for  transverse  component  processing  we  have  a  total  of  nine 
events  to  consider,  for  Rayleigh  wave  processing  this  total  is  ten  events. 

In  ref.  168  transverse  component  processing  we  find  that  the 
waveforms  C  and  D  are  a  common  optimum  for  the  seven  events  167/13, 

184*04,  181*13,  042-05,  047-23,  064-04,  1 54-06;  the  waveform  X  is  the  only 
common  optimum  for  the  events  170/17  and  080-21.  Thus,  optimum  transverse 
component  processing  using  ref.  168  requires  at  least  two  different  reference 
waveform  lengths  for  this  set  of  nine  events.  In  ref.  168  Rayleigh  wave 
processing  the  waveforms  B,  C  an  ^  are  a  common  optimum  for  the  events 
167/13,  170/21,  184*04,  181*13,  064-04;  waveform  A  is  the  common  optimum 
for  170/17,  042-05,  047-23;  Y  is  optimum  for  080-21  and  X  for  154-06.  Thus, 
optimum  Rayleigh  wave  processing  using  ref.  168  requires  four  different 
reference  waveform  lengths  for  this  set  of  ten  events. 

In  ref.  170  transverse  component  processing  waveform 
E  is  the  common  optimum  for  the  eight  events  167/13,  168/15,  184*04,  181*13, 
042-05,  047-23,  064-04,  1 54-06;  X  is  the  optimum  waveform  for  080 -21 . 

Thus,  ref.  170  transverse  component  optimum  processing  requires  two  different 
reference  waveform  lengths.  In  ref.  170  Rayleigh  wave  processing  the  waveform 
C  is  a  common  optimum  for  the  eight  events  168/15,  170/21,  184*04,  181*13, 
042-05,  047-23,  064-04,  156-04;  waveform  E  is  optimum  for  167/13  and  080-21. 
Another  possible  combination  satisfying  optimum  processing  for  this  entire 


set  of  events  is  the  waveform  pair  Y,  D.  Thus,  ref.  170  Rayleigh  wave  optimum 
processing  requires  two  different  reference  wave  form  lengths.  In  this 
respect,  ref.  170  may  be  preferred  over  ref.  168  as  the  optimum  reference 
event  for  this  set  of  events. 

Turning  now  to  ref.  002,  the  optimum  reference  event  for 
the  five  events  281*09,  221*01,  064-08,  084-08  and  060-20,  we  find  that  wave¬ 
forms  C,  D,  E  are  optimum  for  transverse  component  processing  of  the  events 
221*01,  064-08,  084-08  and  waveforms  Y  and  Z  are  optimum  for  281*09  and 
060-20.  In  Rayleigh  wave  processing,  the  waveform  X  is  the  common  optimum 
for  all  five  events.  Thus  ref.  002  transverse  component  processing  requires 
two  different  lengths,  Rayleigh  wave  optimum  processing  requires  only  one 
reference  waveform  length  for  this  set  of  five  events.  Again  we  have  excluded 
results  of  matching  ref.  00  2  with  itself. 

In  most  cases  optimum  processing  requires  that  at  least  two 
different  reference  waveform  lengths  be  generated  from  each  reference 
event  signal  trace.  The  variation  in  optimum  wavefo rm  lengths  indicates  that 
no  general  statement  about  the  optimum  duration  or  signal  content  of  a 
reference  waveform  can  be  made.  However,  examination  of  the  frequency  of 
occurrence  that  a  certain  type  reference  waveform  is  an  optimum,  indicates 
that  some  types  have  a  higher  probability  of  success  than  do  others.  In 
ref*!  168  transverse  component  processing,  types  C  and  D  are  most  frequently 

optimum,  for  ref.  168  Rayleigh  wave  processing  the  most  frequent  optimum 
waveform  lengths  are  A,  B,  C,  D.  For  ref.  170  transverse  component  pro¬ 
cessing,  lengths  C  and  D  are  best,  for  ref.  170  Rayleigh  wave  processing,  C 
and  Y  give  the  best  results.  For  ref.  002  transverse  component  processing, 

C,  D,  E,  Y,  Z  are  best,  and  for  ref.  002  Rayleigh  wave  processing  C,  X,  Y 
prevail.  This  indicates  that  in  general  the  waveform  containing  most  of  the 
signal  bandwidth  and  little  or  no  multipath  arrivals,  such  as  length  C,  may 
have  the  greatest  probability  of  being  the  optimum  reference  waveform  for  a 
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given  reference  event.  Optimum  processing  over  a  number  of  events,  however, 
requires  several  additional  reference  waveform  lengths  per  reference  event. 
There  does  not  appear  to  be  a  straightforward  w  y  to  anticipate  the  additional 
types  of  reference  waveform  required. 

In  the  above  analysis  waveforms  yielding  a  gain  that  is  within 
2  dB  of  the  maximum  gain  are  called  near-optimum  and  compared  on  an  equal 
basis  to  find  the  optimum  waveform  common  to  a  number  of  events.  Relaxing 
this  criterion,  for  instance  to  3  dB,  will  tend  to  reduce  the  number  of  wave¬ 
forms  required;  restricting  the  criterion  tends  to  increase  this  number. 
Obviously,  there  is  a  trade-off  between  the  desired  level  of  optimization,  and 
the  number  of  reference  waveforms  that  such  optimization  requires. 

D.  COMPARISON  OF  OPTIMUM  AND  ROUTINE  PROCESSING  GAINS 

The  optimum  Sinkiang  Province  processing  gains  obtained  in 
this  study  were  compared  with  reported  routine  processing  gains  (Texas 
Instruments,  1972,  1973).  In  routine  reference  waveform  matched  fil¬ 
tering  normally  a  single  high -magnitude  event  is  chosen  as  a  reference 
event  for  a  relatively  large  region,  and  a  part  of  its  signal  trace  which  by 
visual  examination  is  determined  to  contain  most  of  the  signal  bandwidth 
is  taken  as  the  reference  waveform.  In  routine  chirp  waveform  matched 
filtering  chirp  lengths  are  chosen  based  on  distance,  region  and  previous 
processing  experience.  The  optimum  gains  in  this  study  were  obtained  from 
results  of  matched  filtering  using  a  large  number  of  different-length  reference 
waveforms  generated  from  reference  events  at  different  locations  in  the  area 
concerned. 

Table  IV-2  compares  the  optimum  gains  with  the  ALPA 
reference  matched  filter  routine  processing  gains  for  events  processed  in 
both  manners;  Table  IV-3  shows  the  comparison  with  ALPA  routine  chirp 
waveform  matched  filter  gains.  Table  IV-4  presents  a  comparison  of  the 
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ALPA  ROUTINE  REFERENCE  WAVEFORM  PROCESSING 
VS.  OPTIMUM  PROCESSING  GAINS  (dB) 
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Waveform  Matched  Filtering 


ALPA  ROUTINE  CHIRP  WAVEFORM  PROCESSING 
VS.  OPTIMUM  PROCESSING  GAINS  (dB) 
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CENTRAL  ASIA  EVENTS  AVERAGE  ROUTINE  PROCESSING  GAINS  VERSUS 
SINKIANG  PROVINCE  AVERAGE  OPTIMUM  PROCESSING  GAINS  (dB) 
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Chirp  Waveform  Matched  Filtering 


average  optimum  gain  of  groups  I,  II  and  III  events  (excluding  the  reference 
event  002-10)  with  the  average  ALPA  reference  waveform  matched  filter,  ALPA 
chirp  waveform  matched  filter,  and  NORSAR  chirp  waveform  matched  filter 
processing  gains  for  Central  Asia  events  in  the  Sinkiang  Province  and  neigh¬ 
boring  regions  (Sinkiang,  Tadzhik,  Hindu  Kush,  Tibet,  China,  Kirgiz,  Tsinghai, 
Kashmir-Tib^t  Border,  Afghanistan,  Afghanistan-U.  S.  S.  R.  Border).  NORSAR 
reference  waveform  matched  filtering  results  were  not  available  for  Central 
Asia  events. 

The  optimum  processing  method  shows  up  to  8  dB  gain  im¬ 
provements  over  routine  processing.  These  improvements  are  obtained,  how¬ 
ever,  at  the  expense  of  generating  a  large  number  of  different  reference  wave¬ 
forms.  Routine  optimum  processing,  therefore,  would  require  either  laborious, 
or  highly  automated  processing. 

E.  PROCESSING  GAIN  VERSUS  MAGNITUDE,  SIGNAL- TO -NOISE 

RATIO  AND  DISTANCE 

To  complete  the  optimum  reference  waveform  analysis  we 
like  to  investigate  the  influence  of  magnitude,  signal-to -noise  ratio  and  distance 
on  the  matched  filtering  processing  gain.  These  relations  are  shown  in  Figures 
IV-4  through  IV- 12,  per  component,  for  each  reference  event. 

Figures  IV-4,  IV-5,  IV-6  plot  the  matched  filter  processing 
gain  over  bandpass  filtering  versus  event  magnitude,  Figures  IV-7,  IV-8,  IV-9 
versus  event  bandpass  signal-to -noise  ratio,  and  Figures  IV-10,  IV-11,  IV-12 
versus  distance  between  test  event  and  reference  event.  Encircled  data  points 
denote  processing  gains  from  events  for  which  the  reference  event  concerned 
is  optimum.  For  example,  the  gains  of  most  group  I  and  III  events  are  encircled 
on  the  graphs  pertaining  to  refs.  1  6 8  and  170;  most  group  II  events  are  encircled 
on  the  graphs  for  ref.  002. 
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OPTIMUM  GAIN  VERSUS  BODY -WAVE  MAGNITUDE,  REF.  002 
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OPTIMUM  GAIN  VERSUS  BANDPASS  SNR,  REF.  168 
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OPTIMUM  GAIN  VERSUS  BANDPASS  SNR,  REF.  170 
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OPTIMUM  GAIN  VERSUS  DISTANCE,  REF.  168 
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OPTIMUM  GAIN  VERSUS  DISTANCE,  REF.  170 


The  figures  do  not  indicate  any  clear  relation  between  the 
processing  gain  and  the  test  event  magnitude  and/or  signal-to-noise  ratio, 
although  there  is  a  tendency  in  a  few  cases  for  better  gains  at  the  higher  mb's 
and/or  SNR's.  The  only  clear  relation  is  the  trend  of  decreasing  gain  with 
increasing  distance,  as  expected,  for  refs.  168,  170  (Figures  IV-11,  IV-12), 
and  to  a  lesser  extent  for  ref.  002  (Figure  IV-10). 

F.  VERTICAL  AND  RADIAL  COMPONENT  GAIN  DIFFERENCES 

For  a  given  test  event  one  expects  a  certain  type  reference 
waveform  to  act  in  a  similar  manner  for  the  tv  o  Rayleigh  wave  components. 
Where  this  is  not  the  case  the  LR-V  and  LR-R  components  apparently,  in 
addition  to  differing  90°  in  phase,  have  different  signal  and/or  noise  contents. 
From  Table  IV-1  and  Figures  IV-1,  IV-2,  IV-3  we  see  that  these  differences 
occur  rather  frequently.  For  example,  consider  ref.  168  processing  of  event 
170/17  (ref.  170),  where  the  radial  component  shows  a  3.9  dB  better  corre¬ 
lation  than  the  vertical  component.  Table  A-3  (Appendix  A)  shows  the  bandpass 
vertical  and  radial  signals  to  have  a  peak  amplitude  of  464m/r  and  450  m/x  , 
and  RMS  noise  levels  of  3.  5  m fi  and  4.  6m jx  respectively.  Ref.  168,  waveform 
A  matched  filter  processing,  improves  the  signal  amplitudes  to  786  and  828  m/x  , 
the  matched  filter  output  RMS  noise  levels  are  3.0  and  2.  7  mfi  ,  respectively. 
Thus,  where  the  vertical  and  radial  signal  levels  are  improved  with  factors  of 
1.  67  and  1.  84,  or  4.  4  dB  and  5.  2  dB,  respectively,  a  difference  of  0.7  dB,  the 
respective  noise  levels  are  reduced  with  factors  0.  86  and  0.  59,  or  1.  5  dB  and 
4.7  dB,  a  difference  of  3.2  dB,  leading  to  signal-to-noise  ratio  improvements 
of  6.0  dB  for  the  vertical,  and  9.  9  dB  for  the  radial  component.  Thus,  the 
gain  differences  between  the  two  Rayleigh  wave  components  are  established 
for  only  a  relatively  small  part,  0.  7  dB,  by  differences  in  signal  behavior,  and 
for  the  largest  part,  3.  2  dB,  by  the  differences  in  noise  behavior.  This  is 
plausible  since  in  many  cases  the  vertical  and  radial  noise  components  will  be 
at  least  partially  uncorrelated. 
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Similarly,  for  the  events  170/21  (refs.  168,  170  processing), 
181*13  (ref.  168  processing)  and  047-23  (ref.  168  processing),  the  more  than 
2  dB  differences  between  vertical  and  radial  component  processing  gain  are 
mainly  due  to  noise  differences.  However,  the  events  184*04  (ref.  168  pro¬ 
cessing),  064-04  (ref.  170  processing)  and  154-06  (ref.  170  processing )  show 
that  these  differences  are  mainly  caused  by  signal  differences,  possibly  due 
to  imperfect  signal  rotation,  or  multipath  arrivals  with  different  azimuths. 


SECTION  V 
CONCLUSIONS 


In  the  search  for  an  optimum  type  reference  waveform  in 
matched  filter  processing  of  events  of  one  region,  three  major  events  were 
chosen  from  the  Sinkiang  Province  region  as  reference  events.  From  each 
reference  event  signal  trace,  eight  different  reference  waveform  lengths 
were  chosen.  Each  of  these  twenty-four  waveforms  was  applied  as  a  matched 
filter  in  processing  twenty-one  test  events  from  the  same  area.  The  matched 
filter  processing  gains  over  bandpass  filtering  were  analyzed.  The  analysis 
leads  to  the  following  conclusions: 

e  For  the  optimum  reference  event  the  distance  to  the  events 

to  be  processed  is  the  most  important  factor  in  most  cases. 
However,  optimization  of  all  events,  even  within  a  relatively 
small  region,  requires  more  than  one  reference  event.  The 
location  and  characteristics  of  the  additional  reference  events 
are  difficult  to  predetermine.  Optimization  improves  with 
the  number  of  reference  events  available. 

•  The  reference  waveform  length  that  is  optimum  for  the  majority 

of  events  contains  the  entire  signal  bandwidth  but  excludes 
multiple  arrivals.  However,  optimization  of  all  events,  even 
a  small  number  within  a  relatively  small  area,  requires  addi¬ 
tional  lengths.  The  characteristics  of  these  additional  wave¬ 
forms  are  difficult  to  anticipate.  Optimization  improves  with 
the  number  of  reference  waveforms  available. 


•  By  applying  all  of  the  available  24  different  reference  waveform 
lengths  (3  events,  8  lengths  per  event),  2  to  10  dB  gains  re¬ 
lative  to  bandpass  filtering  were  achieved  for  16  out  of  the  21 
test  events  located  over  an  area  of  about  900  x  450  km.  Four 

of  the  five  low-gain  events  were  located  at  least  420  km  from 
the  nearest  reference  event. 

•  Optimum  Sinkiang  Province  processing  yielded  up  to  8  dB 
improvement  over  to  date  routine  processing  of  Central  Asia 
events,  at  the  expense  of  laborious  processing. 

•  No  relationship  could  be  established  between  the  matched  filter 
processing  gain  and  either  the  event  magnitude  or  its  signal-to- 
noise  ratio.  For  two  of  the  three  reference  events  an  obvious 
trend  of  decreasing  gain  with  increasing  distance  was  observed; 
for  the  third  reference  event  a  weaker  trend  was  evident. 

•  Vertical  component  reference  waveform  matched  filter  pro¬ 
cessing  occasionally  gave  more  than  2  dB  different  gains  than 
radial  component  processing  due  to  signal  or  noise  differences 
between  the  two  components. 

Summarizing,  although  significant  processing  gains  are  obtained 
by  a  limited  number  of  predeterminable  reference  event  waveforms,  full  op¬ 
timization  of  reference  waveform  matched  filtering  for  all  events  can  only  be 
obtained  when  a  large  number  of  different  reference  waveform  lengths  is 
generated  from  a  number  of  reference  events  located  throughout  the  area. 

Such  optimization  yields  significant  processing  gains  over  routine  matched 
filter  and  bandpass  filter  processing,  but  requires  either  laborious,  or  highly 
automated  processing.  The  effects  of  such  processing  on  false  alarm  raves 
and  surface  wave  magnitude  measurements  were  not  investigated  in  this  study, 
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APPENDIX  A 


TABULATION  OF  RESULTS 


To  enable  a  detailed  analysis  of  the  results,  the  various  steps 
of  the  signal-to -noise  ratio  improvement  computation  are  combined  in  the 
Tables  A-l  through  A-21,  one  table  per  event. 


Following  the  name  of  the  test  event,  the  first  line  in  each 
table  heading  shows  the  name  of  the  reference  events,  the  second  line  shows 
the  differences  in  great  circle  distance  (D  )  and  azimuth  (A  )  from  array  to 
epicenter,  and  the  great  circle  distance  (S  km),  between  the  test  event  and 
each  reference  event.  The  third  line  indicates  the  waveform  component 
concerned,  LQ-T,  LR-V,  or  LR-R.  The  left  hand  column  in  each  table 
indicates  the  type  of  filter  applied:  BP  for  bandpass,  A  for  type  A  reference 
waveform  matched  filter,  and  so  on.  The  second  column  indicates  the  type 
of  value  presented: 


SIG  = 


NOI  = 


S/N  = 


maximum  value  of  signal  plus  noise  during  the  expected 
signal  arrival  in  m/i. 


RMS  noise  value  computed  over  a  time  gate  just  prior 
to  the  signal  arrival,  in  m/i. 


SIG/NOI,  the  signal-to -noise  ratio  as  defined  in  Section  III. 
D.  3. 

MF  S/N  ,  ,  , 

— — — ~ — -,  the  matched  filter  processing  gain  over  band- 
BP  S/N 

pass  filtering  as  defined  in  Section  III.  D.  4. 


GN 


'  ir-T’-r 


GDB 


I 


=  20  loho 

The  bandpass 


GN,  the  processing  gain  expressed  in  decibels, 
filter  results  SIG,  NOI  and  S/N  are  repeated 


for  each  reference  event. 
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10.7 

NO  I 

1.9 

3.7 

■n 

k9 

3.  R 

1  .n 

l.R 

0 

S/N 

12.3 

12.9 

Hr] 

5.2 

ll.Q 

5 . 6 

GN 

2.0 

1.2 

199 

0.  R 

1  .  1 

7.9 

GCB 

6.1 

1.9 

ms 

HB 

jHMfi 

-1.9 

mg 

1.1 

-1.7 

SIG 

mxmm 

20.  1 

iBHgg 

30.1 

jBQH*' 

NO  I 

WBSk 

2.1 

■ns 

u 

r .  1 

E 

S/N 

9.7 

H9 

9 . 6 

GN 

2.0 

JHI 

1.5 

gtH 

0.9 

■ 

GCB 

6.1 

■Knl 

3.4 

-1.5 

■1 

flSJTl 

-0.7 

IBB 

SIG 

17.6 

24.4 

15.4 

76.? 

JB? 

14.1 

33.4 

14.7 

1 

Enn 

1.9 

2.3 

1.5 

2.3 

I 

1  .6 

3.2 

■ 

BUI 

9.3 

1  C.  6 

7.9 

16.0 

wU 

9.3 

10.2 

4.4 

K|Jj 

1.5 

■gi 

1.6 

1.3 

1.6 

flfl 

1  .  5 

1.0 

1  .0 

BtuH 

3.7 

4.1 

2.3 

3.  R 

■Rl 

3.7 

-0 . 1 

-0 . 1 

■ 

SIG 

Bill 

n  .3 

BiiBJ 

11.0 

21.7 

11.7 

NO  I 

mm 

1.6 

9 

1  .4 

l  .  7 

1  .  P 

S/N 

11.0 

7.  C 

a.O 

12.7 

4.3 

GN 

l.R 

«Oo 

IB 

1.1 

In 

1  .3 

1  .  3 

1  .0 

■ 

GCB 

5.1 

mm 

mu 

1.2 

BOB 

.2.7 

1  .R 

-9.2 

S  IG 

19.9 

74.  1 

2  C  .  0 

n.5 

37.5 

i°.  1 

10.9 

19.0 

13.n 

NO  I 

1 .8 

2.7 

2.0 

3.0 

2.3 

3.9 

1  .4 

1  .  1 

1  .  7 

Z 

S/N 

11.1 

13.6 

10.1 

6 . 6 

1  r.  V 

4.6 

7.9 

1  1  .7 

7  c 

• 

GN 

l.R 

1.2 

1.6 

1.1 

1  .  6 

1.0 

1  .3 

1 . 1 

1  .  3 

GDB 

5.2 

1.7 

3.  P 

0.7 

3.6 

O.  2 

2.2 

1  .0 

1  .  3 

*  C=DEL  T  A  HIFFIOFG)  A=A7IMUTH  OIFF(nFGI  s  =  r>|STANrF  T'~  prpr^v) 


A-  n 


FI  RFF 


S  IN-002-  10AL 


n I F I O-  0.0  A=  O.C  s= 


SIN/16R/1  5AL 


=  42  P 


SIN/ 170/ 1 7  AL 


S  IP, 

314.2 

290.2 

229.6 

217.2 

NO  I 

4.5 

5.P 

4.4 

6.1 

S/N 

69.2 

.  52. P 

52.  P 

35.5 

gn 

3.  1 

1.7 

1.7 

1.6 

GPR 

9.9 

4.4 

4.4 

4.1 

*  F  =  OFL  T  A  niFFlOFGJ  A  =AZ  I  MIJTH  DIFMDFGI  S=niSTANPF  TP  RFF(KR) 


T A R L r  A-  Q 


PROCESSING  RFSILTS,  FVFM  STN-077-07«L 


FI  RFF 
LT  CIF 

fr  ~rrrff 

S  IG 
RP  NO  I 
S/N 
S  IG 
NO  I 

a  s/n 

GN 


S  IN-00 2-  1  Oft L _ 

n=  0.8  A=  1.7  S  =  1C? 

LQ-T  LR-V  LR-P 
24.  8  25. °  22. 5 

9. 4  8.7  8.0 

2.6  8.9  2.° 

38.5  30.  7  3  5.7 

11.8  9.c  1  1  .q 

2.8 
1.1 


SIN/1  6*3/1  5  ft  L 


SIN/! 70/ 1 7  ft | 


*=-1.0  A  =-? .  Q  S  =  ?3ft  |o  =  -Q.q  ft  =  -  p  . y 


75.  7 
11.9 
3.0 
1.1 
C.  5 


LQ-T 
2b,  « 
9.  ft 

3.6 
2  9 . 5 
11.  K 

2.6 


l  R-V 
38.  9 
6.  7 
3.9 

3?.  n 
6.5 
3. ft 
c 


L  P-R 

~ 

3.  0 
3.8 
3  7.5 

10.  1 

3.3 

P 


LO-T 

3ft  .  3 
9.4 
?  •  6 
7  1.6 

11.7 
2.  7 
1  .0 
0.1 


I  P-V 

38.9 
6.7 
3.9 
2  9  .4 
7.4 

4.0 
1 .0 
0.3 


|  P-ft 
33  .q 

3.0 

9.8 

93.0 

1.6 

7  .  1 

l  l 


♦  C  =  r  FI  T  A  r  I  F  F ( n  FG  )  A  =  A7fMUTH  niFF(DFG)  S=HI ST  ANCE  T  D  RFF(KR) 


r 


A-12 


TARLf  A- 1  1 


PROCESSING  «FSLLTS,  FVFNT  STN-0P4-0«QC 


FI 

RFF 

SIN- 

002-1  PM 

MN/1  68/1 5AL 

S15/170/17AI 

LT 

CIF 

n=-i.R 

A=- 1.7 

S  =  268 

0  =—  3 . 6 

A=-c.  4 

S=  6«7 

n=-3 . 5 

A  =  — 5 .6 

5=  6  6  4 

FR 

COM 

Lf.-T 

LR-V 

LR-p 

L0-T 

LR-V 

LP-P 

LC-T 

LP-V 

l_  P  -  W 

SIG 

48. B 

94.2 

77.  6 

4  R  .  0 

94.7 

77.6 

4  P  .  H 

°4  .? 

77.4 

BP 

NO  I 

2.5 

2.4 

2.6 

7.  4 

2.4 

7.6 

2.5 

2.4 

7.6 

■ 

S/N 

19.6 

1  9.6 

30.0 

19.6 

SIG 

gt§,^a 

:iSIK| 

magma 

■ga  9 

79.7 

I 

H 

■  .  1 

wQ 

ns 

tegO 

: 

3.6 

■  '  9 

■  ■ 

v'7; 

0 . 3 

1 

■ 

\ 

■QQ 

9 

0.7 

■ 

■mTI 

■Wfl 

lll§ISicl 

HHf 

-11.1 

SIG 

61.3 

wmm 

60.6 

50.0 

■Ml 

7.5 

Hfl 

7.4 

3.7 

7.6 

B 

17.5 

1  n.  0 

70.7 

77.fi 

mtm 

BOB 

C.Q 

I 

0.9 

0.5 

0.9 

IQ 

Eaa 

■u 

-1  .0 

■MM 

Hi 

-0.7 

-6.5 

-3.4 

SIG 

1 12.7 

154.2 

no.  1 

61. 8 

nn 

■KEgllfj 

i7,Q 

1 

NO  I 

2.6 

1.9 

2.2 

3.3 

7.0 

S/N 

43.0 

80.0 

c  F  .  9 

18. 5 

45.4 

■ 

77.4 

GN 

2.2 

2.0 

2.0 

C.G 

1.7 

0  .« 

■ 

GTB 

6.2 

5.9 

-0.5 

1.3 

mm 

-7.1 

SIG 

IHB 

aEiwni 

97.  7 

94.5 

mmmm 

84 .7 

4  4 . 9 

mm 

■ 

mm 

MB 

7.  1 

7.3 

7.2 

7.7 

1 . 0 

D 

wm 

1 

ii«^i 

43.7 

41.0 

1  9.4 

77.0 

74.7 

■0 

SB 

-  9 

MB 

1.1 

1.4 

0.9 

0.9 

on 

Bmd 

■99 

3H 

C.  9 

7.7 

-0.5 

-0.5 

-l  .7 

S  IG 

iuhi 

■ugj 

& 

«Pgj 

^6.1 

61.1 

51.fi 

NO  I 

mm 

■EH 

*>  r, 

.  • 

7.6 

3.1 

F 

S/N 

■ 

9Q 

\4,’  ■ 

7  4 .6 

GN 

mi 

HS1 

^■ISS: 

s 

BBS 

9 

0  .  P 

wM 

WBB!1 

■99 

mmmsa 

Kp 

-1.7 

m 

SIG 

86.0 

180.8 

144.  C 

71.6 

128.4 

■ftfWl 

46.4 

■ 

2.7 

2.3 

2.3 

2.5 

2.2 

3.0 

7.6 

7  .4 

1.9 

B 

32.3 

g’.-1 

6^.0 

74.4 

5  8.7 

31.9 

7  4.4 

?4 .0 

?  9 . 4 

■ 

1.6 

2.  i 

2.1 

1  .2 

1.9 

1.1 

1  .  ^ 

7. 7 

1.9 

H 

■WWJ 

4.4 

6.6 

6.4 

1.9 

3.5 

0.9 

7.0 

-3 . 7 

-0 . 7 

■ 

SIG 

■J||| 

■ 

■ 

uma 

117. 7 

57.6 

52 .0 

45.9 

NO  I 

mai 

1 

3.7 

7.4 

7.1 

1.4 

S/N 

■ 

36.1 

74 . 1 

79  .  ? 

71.4 

H 

GN 

M 

[Mrs 

1.7 

1  .2 

0.6 

1  .  1 

■ 

GCB 

■99 

mB 11 

mm 

1  .6 

1  .  p 

-3.fi 

o.c 

■ 

S  IG 

89.  4 

157. c 

119.  1 

Pc.7 

125.0 

96.8 

5P  .  P 

76.1 

4  6.fi 

NO  I 

2.4 

2.2 

2.1 

2.4 

2.7 

7.6 

7  .0 

1.6 

B 

S/N 

36.7 

69.7 

c6  •  7 

27.  7 

53.  0 

36.1 

77.2 

?fi  .  ? 

7  1.fi 

GN 

1.9 

l.R 

1  .o 

1.4 

1.4 

1  .  7 

1  .  1 

1  .0 

1 .0 

B 

COB 

5.4 

5  .  C 

c.  5 

7  .  R 

7.7 

1.6 

1  .  1 

-0.7 

0.7 

*  C=PFLTA  DIFFIOEG)  A=A?IMUTH  ^nFintn)  c  =  PTSTANGE  tG  R  F  F  (  K w ) 


Fl|  REF 
IT 

ER|  COM 


SIG 
NO  I 
S/N 


AR L  F  A- 1  2 


PROCESSING  ^F^l.LTSt  EVENT  SI  N*lfl  1  *1  3 /»L 


S IN-002-  1DAL 


"I  STN/1  68/1  SAL 


SIN/170/ 17  A! 


=  2.4  fi= -0 . 7  S 


LO-T  |  LR-V  T  L  P -P 


2.9 

11.0 


GN 

GDR 


S  IG 
NO  I 
R  S/N 
GN 
GOR 


71.8 
ft.  0 
5.7 
1.0 


11.  7 
1.0 

0.2 


GN 

GCR 


SIG 
MO  I 
l  S/N 
GN 
GOB 


22.0 

2.7 

8.2 

1.5 

3.5 


-0.4 


3 


25.6 

2.0 

13.  1 

1.2 

1.5 


2ft. 3 
2.0 

13. C 
1.? 
1.4 


ft.  7 
1.2 
1.3 


15.2 
7.6 
7.  3 
1  .3 

2.  1 


2  F .  8 
7.  1 
<5.7 
1.8 
4.  q 


45.  3 
2.6 
1  7.4 
1.6 

3.  9 


2  5.3 
2.3 
10. 8 
2.0 
5.8 


35  .0 
2.4 
14.3 
1.3 
2.2 


*  r  =  DEL  Ta  OIFF(OEG)  A=A71mUTH  niFF(DFG)  S  =  OISTANCE  TO  R  r  F ( KM) 


?9.0 

7.8 

11.5 

2.0 

ft.O 


TABLE  A-n 


PROCESSING  RFCILTS,  FVFN'T  SU-042-OFQr 


*  C=0ELT A  CIFF(OEC)  A  =A  7  I MUTH  D I  F  F (DFG)  S=niSTANCt  TP  RFF(KW) 


TAB  L  F  A- 1  4 


PROFESSING  R <=  CLI .TS,  FVFNT  S!  N-047-230C 


*  C  =  OEL  T  A  niFF(DFG)  A=A?!MUTH  D I FF  (DFG  )  S=PISTANCF  T  P  R  F  F  (  K  ^  ) 


A-  1  6 


TARLF 


A-l  5 


PPOCFSSING  PFSUTS  FVFNT  5IA-044-04AL 


FI 

REF 

SIN- 

002-10AL 

SIN/1  66/1  5AL 

SI N/l 70/1 7A1 

LT 

fc!F 

0=  2.9 

A  =  3.6 

S  =  494 

9=1.1 

A  =-0.  1 

S=  127 

0=  1.3 

A  =  -0 .2 

s  =  w 

ER 

COM 

LO-T 

LP-V 

LR-R 

LQ-T 

LP-V 

LP-P 

LQ-T 

1  p-v 

( 

StG 

36.  7 

54.6 

44.2 

36.7 

54.6 

44.? 

36.7 

54  .  <• 

44.0 

BP 

NO  I 

4.1 

3.0 

3.5 

4.  1 

3.0 

7.9 

4.1 

3.0 

*>  c 

• 

S/N 

9.0 

1R.2 

12.7 

9  .  C 

1  a.  2 

12.7 

9.0 

1  p.2 

17.7 

TTC 

54 . 4 

R3.1 

80. 0 

41.1 

7  H  .  3 

57.7 

52.1 

72.n 

7  3.9 

NO  l 

4  .  S 

3.5 

4.P 

3.P 

2.  7 

3.3 

4.1 

2.7 

7  .  f: 

A 

S/N 

12.2 

23.6 

16.9 

17. P 

7  R  .  9 

15.7 

12.9 

14.5 

1.4 

GN 

1.3 

1.3 

1.3 

1.4 

1.6 

1.2 

1  .4 

0.8 

0 . 7 

GCR 

2.6 

2.2 

?.  5 

2.7 

4.0 

1.6 

2.9 

-2.0 

-7.6 

SIG 

44.3 

69.7 

69.9 

91.7 

aa.  i 

69.7 

97.2 

91.7 

66 .6 

NOT 

4.6 

3.2 

4.6 

5.8 

7.9 

3.4 

5.3 

7.4 

7.9 

R 

S/N 

9.6 

21.5 

14.6 

15.7 

30.  6 

17.? 

16.0 

.  9 

17.1 

GN 

1.1 

1.2 

1.1 

1.7 

1  .  7 

1.3 

1  .8 

1  .5 

1  .  3 

GOB 

0.5 

1.4 

1.  ? 

4.8 

4.5 

2.6 

5.0 

I  .4 

7.6 

SIG 

55.  * 

49. P 

P9.  6 

57.8 

aa.  5 

69.7 

94.6 

97.4 

5  3.7 

NO  I 

4.2 

3.9 

3.6 

4.9 

3.  C 

3.3 

4.b 

2.3 

7.0 

C 

S/N 

13.2 

14.3 

11.0 

19.9 

2  9.3 

13.2 

19.7 

74 . 6 

I7.', 

GN 

1.5 

0.8 

0.9 

2. 2 

1  .  6 

1.4 

2.7 

1  .9 

1.4 

GOB 

3.3 

-2.1 

-1.3 

6.8 

4.  1 

7.1 

4.7 

c  .6 

7.8 

SIG 

55.1 

57.  3 

50.0 

5  «.6 

84.2 

75.6 

104.9 

7P  .0 

47.7 

NO  I 

4.3 

•>.7 

3.8 

4.6 

3.  1 

7.8 

5.0 

7  .  1 

7.3 

D 

S/N 

12.7 

15.6 

13.  3 

21.7 

76.9 

19.7 

70.3 

75  .  1 

14.9 

GN 

1  .4 

C.5 

l.C 

7.4 

1.5 

1.6 

7.3 

1  .4 

1.7 

GCB 

3.0 

-  1.4 

0.4 

7.4 

3.4 

3 . 0 

7.7 

2  .  p 

i  c 

f  • 

SIG 

53.  1 

57.  1 

46.2 

94.0 

8  5.6 

77. o 

]  07 .6 

77  .6 

5  0.7 

NO  I 

4.0 

3.9 

3.5 

4.4 

3.4 

4  .  Q 

5.1 

T  C 

• 

3.1 

E 

S/N 

13.4 

14.6 

13.1 

21 . 7 

24.  9 

17.9 

7  0.2 

71.0 

14.0 

GN 

1.5 

O.P 

1.0 

7.4 

1.4 

1.4 

2.’ 

1  .  7 

1  .  7 

GOB 

3.4 

-1.9 

0.  3 

7.4 

2.  7 

3.0 

7.0 

1  .  2 

7.0 

SIG 

64.7 

53.9 

27.8 

57.  6 

96.  a 

74.6 

57.5 

04,0 

6  1.7 

NO  I 

4.0 

3.9 

3.5 

5.7 

3.  7 

4.  5 

4  .  P 

7  .  1 

7.4 

X 

S/N 

16.2 

n.9 

8.0 

1C. 9 

?6.0 

16.6 

1  1  .  1 

70 .7 

I7.? 

GN 

l.R 

O.F 

0.6 

1.2 

1.4 

1 . 7 

1.2 

1  .7 

1.4 

GOR 

5.1 

-2.4 

-  4  .  C 

1.6 

3.1 

7.  a 

1  .3 

4.6 

~>  (> 

• 

SIG 

55.6 

5  5.4 

34.  C 

67.4 

90.  5 

78.7 

5  6.6 

84.5 

47.1 

NO  I 

4.1 

3. a 

3.5 

7.6 

3.5 

4.7 

7.5 

7.6 

• 

Y 

S/N 

13.4 

14.5 

5.7 

17.4 

76.0 

1  6.  8 

16.4 

37  .  9 

17.7 

GN 

1.5 

0.8 

0.8 

1.9 

1 . 4 

1  .  7 

1.8 

1  .8 

1  .4 

GCR 

3.4 

-2.C 

-2.  3 

5.7 

7.1 

2.4 

5.2 

5  .  1 

7.7 

S  IG 

50. a 

5  3.  1 

38.2 

74.4 

7«.  9 

72.1 

5  0.9 

6  5.1 

4  7 . 4 

NO  I 

3.7 

’.6 

3.1 

3.8 

7.6 

4.0 

3.3 

2  .  p 

7.7 

Z 

S/N 

n. a 

14.6 

12.4 

15.4 

22.  7 

18.1 

1  »  .  4 

7  7 .0 

14.? 

GN 

1.5 

O.P 

1.0 

7.1 

1  .  ? 

1  .4 

7.0 

1  .  7 

1.7 

GCR 

3.7 

-1.9 

-0.  2 

6.6 

1.7 

3.1 

6.7 

7 .0 

2.1 

*  C=CELTA  0  I F  F  (  0  E  G  )  A  =  4/  I  M(ITF  HIPFIOHG)  S  =  'MSTANrF  yr  i?rrfl<f,) 


A  -  1  7 


T  ft  H  L  £  A  -  1  7 


PRGCFSSING  °F5U  TS»  F  V  r  ft- T  SI  N-l ‘W.-OG  A| 


*  n=PEL  T  A  CIFF(DFG)  A=A7lMUTH  OIFF(PrG>  c  =  n  J  s  T  A'  C  ^  Tf  c  r  p  (  k  ) 


TAPU-  A-l» 


PPnrFSSING  WFSl.l.TS,  FVOT  M  N/l  70/0«  AL 


SIN'/l  6P/1  C&L 


S  !  N / 170/ 1 7  Al 


6.  c 
7.9 
?  .  A 
C.  9 

-1.1 


9.9 

7.7 

’.7 

7.2 

2 

3.7 

2.7 

2.6 

1  .  7 

1  .0 

0  •  Q 

1  0.  7 
l.A 
3.  I 
1  .  1 
0.7 


.  0 
. } 
.7 
.  ? 


10.9 

7.2 

7.7 

1  .  7 


6.7 

7.6 

2.3 

2.7 

7.0 

2.2 

1  .  1 

O.P 

0.7 

*  or  FIT  A  0  IFF  (OFr, )  A  =  AMMUTH  fMrF(PFn)  c=niSTANCF  TO  FFF(KF) 


A-  2 0 


i 


TARLF  A-20 


TAPLF  A-?l 


PROCESSING  sEclLTSt  EVENT  c  I  N-076-2  1  AL 


FI 

REF 

S  IN- 

002-10AL 

CIN/16Q/1  ML 

I  SIN/170/17A| 

LT 

OIF 

0=  4.2 

A  =  0.5 
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